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Abstract 

We present a general algorithm to achieve local 
operators which can produce the GHZ state for an 
arbitrary given three-qubit state. Thus the distil- 
lation process of the state can be realized opti- 
mally. The algorithm is shown to be sufficient for 
the three-qubit state on account of the fact that 
any state for which this distillation algorithm is 
invalid cannot be distilled to the GHZ state by 
any local actions. Moreover, an analytical result 
of distillation operations is achieved for the gen- 
eral state of three qubits. 

PACS numbers: 03.67.-a, 03.65.Bz, 03.67.Hk 

Entanglement manipulation is an important is- 
sue in the studies of quantum information theory. 
On the one hand, it is related to the basic problem 
of which tasks one can accomplish with a given re- 
source of entanglement On the other hand, 
since most applications of quantum information 
theory require the maximally entangled state for 
faithfully transmission of quantum data, it is nec- 
essary to develop the special technique of entan- 
glement manipulation which uses local quantum 
operations and classical communication to purify 
impure entangled states J|, ||. 

There has been extensive work on entanglement 
manipulation. For the case of two-qubit systems, 
the Procrustean method Q provides local oper- 
ations to obtain the maximally entangled state 
from a partly entangled pure state. The gen- 
eral theory of entanglement transformations for 
pure states of the bipartite system has also been 
proposed Q-(7). As far as tripartite states ||- 
[jllj are concerned, researches have shown that 
the GHZ state Jl^] is the maximally entangled 
state which violates Bell inequalities maximally 
and maximize the mutual information of local 
measurements p3| . Hence, it is desirable to pro- 
pose a general algorithm to carry out the local 
operations to distill the GHZ state from a given 
tripartite state. In this paper we present such an 
algorithm. Using this algorithm, we obtain an an- 



alytical result of the distillation operations for the 
general state of three qubits. As far as we know, 
it is the first time for the result of this type to be 
achieved. 

What we address is to transform a tripartite 
state \iPabc) °f three qubits A, B and C into the 
state \i) G Hz) = 75 (I 000 ) + I 111 )) b y local actions 
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where N = (^abc\T{T a ^tIt b ^ c T c \^abc) is 
the probability that the manipulation succeeded. 
Once the operators Ta, Tb and Tq are worked out, 
one can construct the simplest local operations — 
generalized measurements or actions of local fil- 
ters on qubits A, B, and C — to implement the 
distillation process optimally. 

We first give a special representation for the 
three-qubit state \ipABc), which will be useful 
in our distillation method. The representation, 
which we call 'Wootters' representation" here, is 
introduced according to the following idea. Sup- 
pose, the tripartite state \ipABc) considered has a 
form 



ABC) 



(2) 



where {|0a), | 1a)} are the standard basis of qubit 
A, and \^pf~)s are states of the ensemble which re- 
alizes the reduced density matrix p B c of qubits B 
and C. Note that \tpf) is subnormalized, namely, 
the squared length (tpf\ipf) is equal to the prob- 
ability of \ipf) in the ensemble. We then can get 
a kind of representations of \4>abc) by different 
choices of the basis of qubit A 



i=0 



(3) 



with 
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Here i = 0, 1 and the transposition is taken in the 
standard basis {|«a)} (it is included in order to 



be consistent with previous papers, e.g., see Ref. 
PI ). The matrix Va is a representation of the 
unitary operator V A in the standard basis. In- 
troduce the time reversal operation for two-qubit 
states with \(j>) = tr 2 ® tr 2 1 ) , and define the sym- 
metric matrix t,i for {|<^)} with 



m _ 



i = 0,1. 



(5) 



Then according to the studies of Wootters Q, 
there exists a representation 



ABC I 



1 

E 

i=0 



A, 



with 



l 

\A i ) = Ul\iA) = J2(U A ) i j\j A ) 



3=0 



1 

E 



(6) 
(7) 



with 7rf = cannot be distilled to the GHZ state 
by any local actions and classical communication 
acting individually on them. Further discussions 
shall be presented later in this paper. 

Direct observation can be found that the ac- 
tion f R on the qubit R (R = A, B, C) causes the 
reduced matrix of the two other qubits to be sepa- 
rable. For example, consider the operation f a on 
the qubit A. According to the Wootters' represen- 
tation ([j]), after the operation, one shall acquire a 
state (note that such a state is obtained probal- 
istically) with ttq = n^. This means that the 
reduced matrix psc has a zero concurrence, thus 
is nonentangled. 

Theorem 1: A tripartite entangled state of 
three qubits \ipABc) which has Wootters' repre- 
sentations with 7rf ^ (R = A, B, C) can be 
jS^ansformed to the generalized GHZ state by the 
focal operations 



such that 



/a ® /b ® fc\4>ABc) 



(12) 
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Here the parameters irfs, which satisfy ttq > 
> 0, are the square roots of the eigenvalues of 
the Hermitian matrix t^at^ a fl5|| . The difference 
of them defines a "concurrence" fl4| , |l6|| which pro- 
vides a measure of entanglement for the two-qubit 
state pbc- 

Similarly, we can write down two other Woot- 
ters' representations of the state \4>abc) 



where the operators /a, fg, and fc are given by 
Eq. (pi]). Note that a generalized GHZ state 
\ipGGHz) has the following form 
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with relations 
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(9) 



with the parameters < a < (3 by a suitable 
choice of the local basis. 

Proof: Consider the reduced density matrix 
Pbc of qubits B and C. After the local opera- 
tion Ja, it becomes a separable state. Since the 
local operations fs and fc can not produce any 
entanglement for the nonentangled state of qubits 
B and C, it shall still be a separable state in the 
final outcome of the set of actions ( |l2| ) . Accord- 
ing to the communicative property of the set of 
operators {/a, /s, fc}, the same analysis is ap- 
plicable to the reduced density of qubit pairs A-B 
and A-C. Thus, after the local operations (fl2"|), we 
arrive at a state which is pairwise separable. We 
now need to show that a pairwise separable state 



(U R )ij\jR), |xf) = ^(tffl)iil^, three <l ubits is a generalized GHZ state. 



R = B,C, i = 0,l. (10) 



Our distillation algorithm makes use of a set of 
local operations 

fB = (^) 1/2 \R )(Ro\ + \Ri)(Ri\ 



= U R [(^)y 2 \0 R )(0 R \+\l R )(l R \]U R ,R = A,B,C. 

(11) 

Note that these operations require ir^- 7^ since 
for the case of ir^ = 0, f R becomes a projection 
operator of the pure state of qubit R and it shall 
disentangle the tripartite state. In fact, the states 



Lemma 1: A tripartite entangled state of three 
qubits is pairwise separable if and only if it is a 
generalized GHZ state. 

Proof: It is obvious that a generalized GHZ 
state (|l^) is pairwise separable. To prove the con- 
verse statement, we use the fact that an arbitrary 
tripartite state of three qubits can be written in 
the form || 

\4>abc) = Ao|000)+A 1 e lv |100)+A 2 |101)+A 3 |110)+A 4 |lll) 

(14) 

by a suitable choice of local basis, where the co- 
efficients A; are all real and non-negative and 
tp is a phase between and w. This represen- 
tation is a minimal description, with only five 
terms, for the three-qubit states. Now, sup- 
pose a tripartite entangled state IVabc) w hich 



has a minimal representation with coefficients Af 
(i = 0, ■ ■ ■ , 4) and a phase tp s is pairwise separa- 
ble. Then, due to separability, there are relations 
(R = A.B.C) for the Wootters' repre- 



sentations of \4>abc)- 
tion yields 



A straightforward calcula- 
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Noticing that \iPabc) i s a tripartite entangled 
state, i.e., it can not be written as a bisepara- 
ble form and a product form of the three qubits, 



we then have the solutions A? 



An 



for 



the above set of equations. Now the state IV'abc) 
takes the form 



IV'W) = AglOOO) + AI|111> 



(18) 



This completes our proof of Lemma 1 and then 
Theorem 1. 



Let |V>' 



denote the state of the outcome 



ABC/ 

corresponding to the local operations of (|12|) . Ac- 
cording to Theorem 1, we have 

|VW=«|0W>+/3|1'1'1'). (19) 

The coefficients a, j3 and the local basis {|0') ; |1')} 
of the three qubits are completely determined by 
the initially given state \tpABc}- They can be di- 
rectly calculated from Eqs. (|ll|) and (|l2|). Now 
to acquire the GHZ state one only needs to carry 
out one of the following actions (noticing that we 
have set a < (3) 



l' R ){l' R \, R = A,B,C (20) 



on the qubits A, B or C respectively. 

The algorithm established above shall enable 
one to achieve local operators of the distilla- 
tion actions for any pure tripartite state of three 
qubits. Before presenting the general result for 
such distillation actions, we give some discussions 
which our algorithm implies. 

In general, to distill the GHZ state from a 
three-qubit state requires local operations on each 
of the three qubits. Nevertheless, there are excep- 
tions. The simplest example is the generalized 
GHZ state for which one only needs to perform a 
local operation on any one of the three qubits. In 
detail, on which qubits the operations to be per- 
formed is determined by properties of the three 
reduced matrices pab, Pac, an d Pbc — whether 
they are separable or not. For instance, assume 
Pab of the state \ipABc) is separable. Then we 
have ttq = 7rf due to its separability. This results 
in that the local operator fc of Eq. (nil) becomes 



an identity operator. A noticeable example of this 
type is the slice state |0 which has the form 



A |000) + Ai|100) + A 4 |lll). (21) 



A simple calculation shows that the reduced ma- 
trices pab and pac of it are separable. Thus to 
achieve the GHZ state from \tp s u ce ) one only needs 
to perform a local operation Ta on the qubit A. 

Now we present an analytical result of distilla- 
tion operators Ta^Tb <8 Tc in ([!]) for the gen- 
eral pure three-qubit state. We still use the mini- 
mal representation (|lj) for the state \ipABc)- Di- 
rect calculations give the parameters of the state 
\iPabc) 



< 1 = J\A\* + \l\l±\A\ 



<i =A (VA| + Al±A 2 ), 



'0,1 



(22) 



where '+' is for ttq and '— ' for 7r^. The complex 
number A is defined as 



A = A 2 A 3 - A^e* 



(23) 



The local operations which transform the state 
lipABc) into the GHZ state are given by 

Ta®T b ®T c = U A f A ® U' b Jb ® U'cfcfa (24) 

where the operators are shown as follows: 

f R : f R (R - A,B, C) is given by Eq. ©. The 
unitary operator U R in it comprises two parts 



(25) 



Here, Urq is a unitary transformation which di- 
agonalizes the Hermitian matrix r„j?r*«. Direct 
calculations give that 



U A0 



where 
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a i. ^2 — V A 2 + A4 A3 — >/A§ + A 

Ob = arctan , 0c = arctan ^j— 

A4 

The operator 

UlU=i\l R )(lR\ + \0R)(0R\ 



(27) 



(28) 



is included so that the diagonal elements of 
Urt^rU^ become real and non-negative. 

f' R : The operator f' R (R = A, B, or C) of 
(20) transforms the generalized GHZ state into 
the GHZ state. It is equivalent to perform any 
one of them. In equation (E4h we have chosen the 



operation f' c acting on the qubit C. To achieve 
f' R requires the detailed knowledge of parameters 
and local basis of the generalized GHZ state fll9| ) . 
After some hard but straightforward calculations, 
we obtain the ratio of the two coefficients 



constant under the invertible local operations of 
qubits B and C (see Ref. |ll| the relation 



(A| + Ai)(Al + Ai) 



-,1/2 



|A| 2 + A^AI 
and the local basis of the qubit C 
V2i, 



(29) 



10; 



cl 



-[(cos? 



y/2i„ 
= — [(cos I 



?c)\0c 



'c - sin 9c ) | Oc ) - (cos 9c +sin 9 c ) 1 1 

(30) 
+ AqA|, 
E3I) and 



Here we assume (A 2 + A|) (A| + A|) < 
and the parameter A is given in Eq. 
in (|7|). 

The final local unitary transformations 



U'r- 
U' R 



\QrWr\ + \U)W r \ R = A,B,C (31) 



7rf = Tii — tti — shall be retained under 
the actions Tb ® Tc- Similar analysis gives that 
this relation shall also hold under the action of 
Ta- Thus we can conclude that local operations 
Ta®Tb® Tc shall take a W-class state to another 
W-class state, so that the GHZ state will never be 
produced. 

In summary, we have presented an algorithm 
to distill the GHZ state from a single copy of the 
c)|lctHree-qubit state. It enable one to achieve directly 
the operators of the distillation operations, thus 
the distillation process can be realized optimally. 
We then apply our distillation algorithm to the 
general state of three qubits and obtain an ana- 
lytical result of operations for such a process. Fi- 
nally, we show that the state for which our distil- 
lation algorithm is not valid cannot be distilled to 
the GHZ state by any local operations. Thus the 
distillation algorithm we presented is sufficient for 
the tripartite state of three qubits. 



c 



are included to revert the local basis {|0 fl ), 
of generalized GHZ state ( |l9| ) to the initial 
standard basis {|0h),|1r)} of the state \ipABc)- 



Acknowledgments 



{loy, |l' c )} have been given in Eq. (|30|). Now we 
present the local basis of qubits A and B 



|0a> = ^[(cos0 2 -t-sin0 2 )e- 



ze3 |(U) + (cos( 



This work was supported in part by the Na- 
tional Natural Science Foundation, the Doctoral 
Education Fund of the Education Ministry, and 
the Nuclear Theory Fund of HIRFL of China. 
sin# 2 ) e ^|U)], 



II'.) = 



V2i 



[(cos I 



10' ) = 



V2i.. 
= —[(cos i 



|1 B ) = — [(cos0 B -sin 



References 

2 -sm9 2 )e~ l93 \0 A )~{cos92+sm9 2 y 93 \lA)], 

[1] C.H. Bennett, Phys. Today 48, No. 10, 24 

(1995) . 

B +sin 6» B )|0 s ) + (cos 6» B -sin 9 b )\1b)], 

[2] C.H. Bennett, H. Bernstein, S. Popescu, 
and B. Schumacher, Phys. Rev. A 53, 2046 

(1996) . 

(32) 



>B)|O B )-(cos0 B +sin0 B )|l s > 



where 9 2 , #3, and 9b are given by (|27|). 

There remains one case to consider, namely, the 
states with parameters 7rf = (R — A, B, C) for 
their Wootters' representations. For these states 
our distillation algorithm is not valid any more. In 
fact, as far as the tripartite entangled state is con- 
cerned, the three relations 7r^ — for R — A, B 
and C are equivalent. For example, the condition 
irf = implies that 7rf = 7rf = 0. This can be 
seen from Eq. (p2|). Assume irf = for the state 
\ipABc)- Since Ao ^ (otherwise the state \4>abc) 
will be biseparable), the relation = leads to 



[3] G. Vidal, Journ. of Mod. Opt. 47, 355 (2000) 

[4] 



[5] 



C. H. Bennett, G. Brassard, S. Popescu, B. 
Schumacher, J. A. Smolin, and W.K. Woot- 
ters, Phys. Rev. Lett. 76, 722 (1996). 

M.A. Nielsen, Phys. Rev. Lett. 83, 436 
(1999). 



A4 = 0, and then 7rf 
now has the form 



b - „c _ 



0. The state \^>abc) 



[6] G. Vidal, Phys. Rev. Lett. 83, 1046 (1999). 

[7] D. Jonathan and M.B. Plcnio, Phys. Rev. 
Lett. 83, 1455 (1999). 

[8] O. Cohen and T.A. Brun, Phys. Rev. Lett. 
84, 5908 (2000). 



\^abc) = A |000)+A 1 e^|100)+A 2 |101)+A 3 |110). 

(33) 

This state is the so-called "W-class state" JTof 
which could not be distilled to the GHZ state 
by any local actions. It can be understood in 
the following way. Since, the ratio tti/ttq is a 



[9] 



[10] 



A. Acin, A. Andrianov, L. Costa, E. Jane, J.I. 
Latorre, and R. Tarrach, Phys. Rev. Lett. 85, 
1560 (2000). 



W. Diir, G. Vidal and J.I. Cirac, |quant- 
ph/0005115[ 



[11] T.A. Brun and O. Cohen, |quant-ph/0005124 



[12] D.M. Greengerger, M. Home, A. Zeilinger, 
Bell's theorem, Quantum Theory, and the 
Conceptions of the Universe, ed. M. Kafatos, 
Kluwer, Dordrecht 69 (1989). 

[13] N. Gisin and H. Beschmann-Pasquinucci, 
Phys. Lett. A 246, 1 (1998). 

[14] W.K. Wootters, Phys. Rev. Lett. 80, 2245 
(1998). 

[15] R.A. Horn and C.R. Johnson, Matrix Analy- 
sis (Cambridge University Press, New York, 
1985), p. 205. 

[16] S. Hill and W.K. Wootters, Phys. Rev. Lett. 
78, 5022 (1997). 

[17] H.A. Carteret and A. Sudbery, J. Phys. A 33, 
4981 (2000). 

[18] N. Linden, S. Massar, and S. Popescu, Phys. 
Rev. Lett. 81, 3279 (1998). 

[19] A. Kent, N. Linden, and S. Massar, Phys. 
Rev. Lett. 83, 2656 (1999). 



